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Magnetization steps (MST's) from Mn 2+ pairs in several single crystals of Zni-^Mn^O 
(0.0056<x<0.030), and in one powder (x=0.029), were observed. They were used to determine 
the four largest exchange constants (largest J's), and the single-ion axial anisotropy parameter, D. 
The largest two exchange constants, Ji/fcB=— 18.2±0.5 K and J[/ftB~~ 24.3±0.6 K, were obtained 
from large peaks in the differential susceptibility, dM/dH, measured in pulsed magnetic fields, H, 
up to 500 kOe. These two largest J's are associated with the two inequivalent classes of nearest 
neighbors (NN's) in the wurtzite structure. The 29% difference between Ji and J[ is substantially 
larger than f3% in Cdi-iMnxS, and 15% in Cdi-^Mn^Se. The pulsed-field data also indicate that, 
despite the direct contact between the samples and a superfluid-helium bath, substantial departures 
from thermal equilibrium occurred during the 7.4 ms pulse. The third- and fourth-largest J's were 
determined from the magnetization AI at 20 mK, measured in dc magnetic fields H up to 90 kOe. 
Both field orientations H || c and H || [1010] were studied. (The [lOlO] direction is perpendicular to 
the c-axis, [0001].) By definition, neighbors which are not NN's are distant neighbors (DN's). The 
largest DN exchange constant (third-largest overall), has the value J/fca=— 0.543±0.005 K, and is 
associated with the DN at r=c. Because this is not the closest DN, this result implies that the 
J's do not decrease monotonically with the distance r. The second-largest DN exchange constant 
(fourth-largest overall), has the value J/fcB~— 0.080 K. It is associated with one of the two classes 
of neighbors that have a coordination number z n — 12, but the evidence is insufficient for a definite 
unique choice. The dependence of M on the direction of H gives D/kB~— 0.039±0.008 K, in fair 
agreement with —0.031 K from earlier EPR work. 

PACS numbers: 75.50.Ee, 71.70.Gm, 75.10.Jm, 75.60.Ej 



I. INTRODUCTION 

The most extensively studied diluted magnetic semi- 
conductors (DMS's) are IF VI materials (A n B VI , where 
A=Zn, Cd, B=S, Se, Te) in which some of the cations 
have been replaced by manganesei The magnetization- 
step (MST) method is one of the most effective tech- 
niques of measuring antiferromagnetic (AF) exchange 
constants in DMS'si^ This technique has been used to 
determine nearest-neighbor (NN) and distant neighbor 
(DN) exchange constants in several II- VI DMS's with 
the zinc-blend o 4 i 5 ' 6 and wurtzit o 7 ! 8 ' 9 structures. Rele- 
vant theoretical treat ments of these exchange constants 
include those in Refs. lull lll2l In addition to exchange 
constants, the MST method gives information about 
magnetic anisotropies, and about the distribution of the 



magnetic ions in the crystal, on a length scale of several 
atomic dimensions. 



A new class of II- VI DMS's based on ZnO, especially 
Ziii-^Mn^O, has attracted attention recently because 
theoretical calculations suggested the possibility of fer- 
romagnetism above 300 K in p-type samples^ Experi- 
mental works on epitaxial thin films of Zni-^Mn^O gave 
different results: ferromagnetism was reported in Ref.ll4l 
but according to Ref. |lja the largest exchange constant 
is antiferromagnetic, J/fcB~— 15 K. In the present work, 
MST's from several Zni-^Mn^O single crystals, and from 
one powder sample, were used to determine the largest 
four exchange constants. The single-ion axial anisotropy 
parameter D was also determined. 
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FIG. 1: The wurtzite crystal structure. The large white 
spheres are the cations, the small black spheres are the anions. 
The "central cation" is labeled as X. One example of each of 
the neighbor classes in Table [I] is indicated by the number n 
specifying that class. 



TABLE I: Classification of neighbors in the vicinity of a "cen- 
tral cation" in the hep cation structure. The neighbor class 
is specified by n. The distance of such a neighbor from the 
central cation, in the ideal hep structure, is r n . The coordi- 
nation number z n is the number of neighbors of class n which 
surround the central cation. The fourth row gives alterna- 
tive designations for the exchange constants J(n), e.g., the 
exchange constant J(4) for a neighbor of the class n=4 is des- 
ignated as J3. The superscripts "in" and "out" distinguish 
between equidistant but inequivalent neighbors: those in the 
same c-plane and those in different c-planes. The dipole- 
dipole interaction constant (gus) 2 /r^, expressed in kelvin, is 
for the lattice parameter a of ZnO, but using the ideal ratio 
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II. CRYSTAL STRUCTURE AND 
CLASSIFICATION OF NEIGHBORS 

A. Classification of Neighbors by Classes 

The hexagonal wurtzite structure of ZnO (space group 
P63TOC) is shown in Fig. ^ The cations, (open circles) 
form an hexagonal close-packed (hep) structure. The 
cation marked by X is chosen as the "central cation." The 
other cations are often classified by their distances from 
the central cation, i.e., nearest neighbors (NN's), second- 
neighbors, etc^fi A major shortcoming of the classifica- 
tion by distance, is that in the (ideal) hep structure some 
equidistant cations are not equivalent from symmetry 



point of view. Equidistant but symmetry-inequivalent 
cations have different isotropic exchange constants. 

The classification of neighbors by symmetry, instead 
of distance, is discussed in Refs. andU In this clas- 
sification, neighbors are divided into "classes." Neigh- 
bors of the same class have the following property: When 
cation sites are occupied by magnetic ions, neighbors of 
the same class have the same isotropic (Heisenberg) ex- 
change interaction with the magnetic ion at the central 
site. The exchange constant J is therefore the same for 
all neighbors of the same class. The underlying reason is 
that all pairs of cation sites consisting of the central site 
and a neighbor of a given class are related to each other 
by operations of the space group of the cation structure. 
Interactions other than isotropic exchange sometimes re- 
quire distinctions between neighbors of the same class^ 

Properties of several classes of neighbors are given in 
Table [I] The number n, which is the same as in Ref.H 
specifies the neighbor class. Note that n=l and n=2 
are two inequivalent classes of NN's. As can be seen in 
Fig. ^ n=l corresponds to "in plane" NN's (i.e., NN's 
which are in the same c-plane), whereas n=2 corresponds 
to "out of plane" NN's. The distances in Table [I] are for 
the ideal hep structure, with c/a—y/8/3, where a is the 
NN distance. A neighbor of the symmetry class n=3 
corresponds to a second neighbor in the classification by 
distance. Neighbors of class n=4 are reached from the 
central cation by moving a distance c along the c-axis. 
They are the closest neighbors along the hexagonal di- 
rection. Neighbors of classes n=5 and n=6 are equidis- 
tant but are inequivalent by symmetry. The remaining 
neighbor classes in Table [I] n~7 and n=8, are included 
in Fig. ^ The parameter z n in Table [I] is the "coordi- 
nation number," i.e., the number of neighbors of class n 
surrounding the central cation. 

The notation for the exchange constants J's associated 
with different neighbors has been evolving, to accommo- 
date newer classifications of these J's. In early works, 
when neighbors were classified by their distance r, the 
notation for the J's also was based on distance: Ji for 
NN's, J2 for next-nearest(second)neighbors, J3, for third 
neighbors, etc. For II- VI DMS's with the zinc-blende 
structure (fee cation lattice) this notation is still quite 
useful because each of the eight shortest distances r is 
associated with a unique neighbor class n (see footnote 
115 in Ref. 0). However, the distance-based notation is 
totally inadequate for DMS's with the wurtzite structure 
(hep cation structure). In the ideal hep structure, the 
shortest distance r already corresponds to two classes of 
NN's, with different J's. 

An early apparent advantage of the distance-based no- 
tation followed from the prediction^ that the magnitudes 
(sizes) of the J's decrease monotonically with increasing 
r. If true, this prediction would have made the classifi- 
cation by distance equivalent to a classification by size. 
However, later theoriesjii*i£ and recent experiments;**^ 
indicate that there is no simple correspondence between 
size and distance. 
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The different notations for the J's that are used in the 
present work serve different needs. The simplest nota- 
tion, J(n), associates the J's with the neighbor classes n 
listed in Tabled For example, J(4) is the exchange con- 
stant with a neighbor of the class n=4. The disadvan- 
tages of this notation is that neither the relevant distance 
r nor the ranking by size are immediately obvious. 

An alternative notation, similar to that in Ref. Il8l is 
given in the fourth row of Table [I] This notation too 
is based on division of neighbors into symmetry classes, 
but it also gives some information about the distance r. 
Instead of the number n, the neighbor class is specified 
by a combination of subscripts and superscripts. The in- 
formation about distance is given by the numerical value 
of the subscript, which increases with increasing r. For 
the same subscript, the superscripts "in" and "out" are 
used to distinguish between equidistant but inequivalent 
neighbor classes. Thus, JJ n and J° ut are the J's for the 
two classes of NN's. A prime is added as a superscript 
to indicate that the distance r is approximately, but not 
exactly, the same as for an unprimcd exchange constant 
with the same subscript. For example, the exchange con- 
stants for the neighbor classes n=7 and n=8, whose dis- 
tances differ by only 4%, are designated as J' A and J4, 
respectively. 

B. Classification of Exchange Constants by Size 

The preceding two notations for the J's were both 
based on the division of neighbors into symmetry classes. 
This classification, however, has a serious practical draw- 
back. Quite often the magnitude (size) of an exchange 
constant is measured before the neighbor class with which 
it is associated is determined. Prior to such a determi- 
nation, any notation based on the neighbor class is not 
useful. It is then more practical to adopt a notation that 
is based primarily on the ranking of the J's by size. In 
the present work, only the four largest J's were mea- 
sured. The chosen designations of these J's in terms of 
their sizes are as follows: 

1. The largest two exchange constants are labeled 
Ji and J{, with J\ chosen (arbitrarily) to be the 
smaller of the two. These two exchange constants 
are associated with the two inequivalent classes of 
NN's. In the present work it has not been deter- 
mined which of the two corresponds to JJ n and 
which to Jf ut . 

2. By definition, any exchange constant J which is not 
associated with either of the two classes of NN's is 
a distant-neighbor (DN) exchange constant. The 
largest DN exchange constant (third-largest over- 
all) is called j( 2 >, the second-largest DN exchange 
constant (fourth-largest overall) is called j( 3 \ 

The assignment of and to specific neighbor 

classes is a major task class that will be discussed in 
detail. 



III. EXPERIMENTAL 
A. Samples 

Single crystals of Zn^.Mn.O (0.0056<x<0.030) were 
grown by chemical vapor transport using chlorine as the 
transporting agent. The growth temperature was 900°C. 
The Mn concentration x was obtained using three meth- 
ods: 

1. From the Curie constant, obtained from a fit of 
the susceptibility between 200 and 300 K to a sum 
of a Curie- Weiss susceptibility and a constant \d 
representing the lattice diamagnetism. The values 
5=5/2 and 5=2.0016 (Ref. 19) for the Mn 2 + ion 
were used. 

2. From the apparent saturation value M s of the mag- 
netization. The determination of M s (also known 
as the "technical saturation value" ) is discussed in 
Sec. [V] The relation between x and M s was dis- 
cussed in Refs. I2I andl2pl among others. This rela- 
tion is based on the assumption of a random distri- 
bution of the Mn ions over the cation sites. 

3. From atomic emission spectroscopy with induc- 
tively coupled plasma (ICP-AES). 

Table |n] compares the results of the three methods. 
The good agreement indicates that the apparent satura- 
tion value M s is consistent with a random distribution of 
the Mn ions. The last column in Table ITTI gives the chosen 
values of x that will be used to label the various samples, 
namely, £=0.0056, 0.021, 0.029, and 0.030. These values 
will also be used in the data analysis. 

X-ray powder diffraction data were obtained on small 
portions of the samples with x—0. 0056, 0.021 and 0.030. 
These data were taken with a Bruker model "D8 Discover 
with GADDS" spectrometer, using Cu-K Q radiation. All 
the powder diffraction patterns were in good agreement 
with the wurtzite structure (space group P&^mc). No 
other crystallographic phase was detected. 

The samples used in measurements of the magnetiza- 
tion M, and of the differential susceptibility dM/dH, had 
linear dimensions of 2 to 4 mm. The only exception 
was one set of pulsed-field data on a powder obtained 
by crushing the single crystal with £=0.029. Pulsed-field 
data on that single crystal were obtained before it was 
crushed. 



B. Magnetization measurements 

Three types of magnetization measurements were per- 
formed: 

1. The magnetization M at T=20 mK was measured 
in dc magnetic fields up to 90 kOe. These data 
were taken with a force magnetometer operating 
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TABLE II: Mn concentration, x, as obtained from: 1) the ap- 
parent saturation value M s , 2) the susceptibility between 200 
and 300 K, and 3) atomic emission spectroscopy with induc- 
tively coupled plasma (ICP-AES). The "best value" (B.V.) is 
the value adopted in the text. 



Sample 
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in a plastic dilution refrigerator. The experimental 
techniques were described earlier^i The magnetic 
field H was either parallel to the c-axis, or parallel 
to the [1010] direction (one of the directions perpen- 
dicular to the c-axis.) The dc magnetic field was 
produced by a NbTi superconducting magnet. Sev- 
eral traces of M vs. H, in both increasing and de- 
creasing H, were taken for each experimental con- 
figuration. All such traces were similar, and showed 
no hysteresis. They were averaged in order to im- 
prove the signal/noise ratio in the final result. 

2. The magnetization M was measured at 0.65 K in dc 
fields up to 170 kOe. The samples were immersed 
in a liquid 3 Hc bath. A vibrating sample magne- 
tometer (VSM) operating in an 18-tesla supercon- 
ducting magnet (Nb 3 Sn wire) was usedi 2 ^ Again, 
no hysteresis was observed, and traces were aver- 
aged to improve the signal/noise ratio. 

3. The differential susceptibility dM/dH was mea- 
sured in pulsed magnetic fields up to 500 kOe. The 
experimental techniques were described earlier^ 
The pulse duration was 7.4 ms. The sample was 
in direct contact with a liquid 4 He bath main- 
tained at a temperature Tb a th=l-5 K. However, 
the data showed that despite the direct contact, 
the sample was not in thermal equilibrium with 
the liquid-helium bath during the pulse. Such non- 
equilibrium effects were found in earlier pulsed-field 
experiments. 23 ! 24 ' 25 ! 26 



IV. NN EXCHANGE CONSTANTS FROM 
PULSED-FIELD DATA 

A. NN Exchange Constants 

Figure [21 shows dM/dH vs. H for the single crystal 
with 2=0.029. This trace is from the field-down por- 
tion of a pulse with a maximum field of 500 kOe. The 
two prominent peaks at high fields, labeled as 1 and 1', 
correspond to two MST's. The expanded view of these 
two peaks, shown in Fig. |3 indicates that both peaks 
have similar heights and widths. The maxima are at 
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FIG. 2: Differential susceptibility dM/dH for a single crystal 
with £=0.029, measured in pulsed fields. These results are 
for the "down" portion of the pulse (decreasing H), with the 
magnetic field parallel to the c-axis. The two large peaks, 1 
and 1', are attributed to the two inequivalent NN pairs. The 
small peaks Q and F are discussed in the text. 
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FIG. 3: The high field portion of the results in Fig. H The 
raw experimental data, from Fig. [5] are represented by the 
dotted curve. The solid curve was obtained by subtracting 
a linear baseline. The two dominant peaks 1 and 1' are at- 
tributed to the two inequivalent NN pairs in the wurtzite 
structure. 



H i=275 kOe and H[=367 kOe. These two peaks were 
not resolved in the up portion of this field pulse. 

Pulsed-field measurements were also performed on sin- 
gle crystals with £=0.021 and 0.030, and on a powder 
obtained by crushing the single crystal with £=0.029 (af- 
ter the data in Figs. Inland |21 were obtained). In all cases 
the large peaks, 1 and 1 , were resolved in the down por- 
tion of the pulse. However, in the up portion of the pulse 
these peaks were well resolved only in the following situa- 
tions: a) for the powder sample, in all field pulses, and b) 
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for the single crystal with £=0.029 when the maximum 
field was 420 kOe, as compared to 500 kOe for the pulse 
in Fig. |2 These results are explained later. 

Based on the data in all samples, the two large MST's 
observed in pulsed fields are at f/i=270±8 kOe and 
-ffj=362±8 kOe. Because the uncertainties in these two 
values are correlated, the uncertainty in the difference is 
much smaller, i.e., (H[ — ffi)=92±2 kOe. The average is 
Hx={H[ + Hi)/2=316±8 kOe. 

Earlier data for other DMS's with the wurtzite 
structure 7 ^ showed that: a) each MST from NN pairs 
splits into a doublet, corresponding to the two inequiv- 
alent classes of NN's, and b) the two corresponding ex- 
change constants, J[ and J\, are the largest. The two 
large MST's at H[ and Hi are therefore attributed to 
the two inequivalent NN pairs. 27 

For Mn 2+ pairs with intra-pair exchange constant J, 
the magnetic fields H n at the MST's are given by^ 

gf x B H n = 2n\J\, (1) 

where n=l, 2, . . . , 5. In the present case the calculated 
deviations from Eq. (JTJ, caused by anisotropics and DN 
interactions, turn out to be smaller than the experimental 
uncertainties in Hi and H[. Using n=l and c/=2.0016, we 
obtained Ji/fe B =-18.2±0.5 K and J 1 /fc B =-24.3±0.5 K. 
Unfortunately, it was not possible to conclude which of 
the two exchange constants is JJ n and which is J° ut . In 
the case of Cdi-sMii^Se^ the smaller of the two NN 
exchange constants, defined here as Ji, was identified 
as J° ut , and the larger as J[ n . This identification was 
based on the effect of the Dzyaloshinkii-Moriya (DM) in- 
teraction on the widths of the MST's, and it also agreed 
with Larson's prediction^ In the present work the ef- 
fect of the DM interaction was not apparent in the data, 
presumably because this interaction decreases rapidly as 
the atomic number of the anion decreases, i.e., much 
smaller for oxygen than for selenium. 29 Although direct 
evidence is lacking, based on the experimental result for 
Cdi-nMn^Se and on the theory, we speculate that in the 
present material too, Ji—J° ut and J[—J[ n . 

As discussed earlier^ the dominant superexchange 
path for both classes of NN's, which is through the in- 
tervening anion, is the same. The difference between 
JJ n and J° ut is attributed to differences in the other 
exchange paths whose contribution is smaller. In the 
present material AJi=| J[ n ~ J° ut \ is 29% of the aver- 
age Ji/fcs=— 21.2±0.5 K. This percentage difference 
should be compared to 13% in Cdi-^Mn^S, and 15% 
in Cdi^Mn^Sei 7 ^ 



B. Other Features of the Pulsed-Field Data 

Figure OH also shows two small peaks, labeled as Q and 
F. Computer simulations indicate that these peaks are 
due in part to MST's from quartets (tetramers). Each 
of these quartets consists of four spins that are coupled 
by some combination of Ji and/or J[ exchange bonds. 



Other possible contributions to these small peaks may 
be due to cross-relaxation processes that can occur in 
the absence of equilibriumiSSi2L2i The non-equilibrium 
behavior in pulsed fields, including cross- relaxation, will 
be discussed in a later publication. Very briefly, the cross 
relaxation processes that may contribute to Q are simi- 
lar to those discussed in Ref. |3jJ in connection with the 
second- harmonic peak Pi/2- The cross relaxation pro- 
cess that may contribute to F involves two NN pairs 
of different classes. The latter process is most rapid at 
H\={H'i + Hi)/ '2, where the energy separation between 
the two lowest levels for one class of NN pair matches 
that for the other class of NN pair. 

The widths of peaks 1 and 1' provide convincing ev- 
idence for other types of non-equilibrium processes, as- 
sociated with a phonon-bottleneck, which restricts the 
sample-to-bath heat fiowi 24 ! 2 '?!? 2 !??!? 4 ^ From Fig. the 
full width at half height of cither of these two peaks is 
14 kOe. This value should be compared with a minimum 
equilibrium width of 39 kOe at Tb a th=l-5 K. In addition 
to narrow widths, peaks 1 and 1' also show a pronounced 
asymmetry which is characteristic of non-equilibrium be- 
havior resulting from a phonon bottleneck. The 14 kOe 
width at half height is the sum of 5.5 kOe from the rise 
and 8.5 kOe from the fall. (The data in Fig.|3are for de- 
creasing H, so that the rise corresponds to fields above 
the maximum of dM/dH, and the fall corresponds to 
fields below the maximum.) Both the narrowing and the 
asymmetry are predicted from models for the phonon 
bottleneck. 

The phonon bottleneck also accounts for the difficulty 
of resolving the peaks I and 1' during the up portion of 
the field pulse. At the beginning of the pulse the mag- 
netocaloric effect associated with the alignment of the 
singles causes the sample to warm. This heating is basi- 
cally the inverse of cooling by adiabatic demagnetization, 
except that it is not fully adiabatic. To observe the peaks 
I and 1', much of this heat must be transferred to the 
bath before these peaks are reached. In the up portion of 
the pulse the peaks are reached earlier than in the down 
portion. Apparently, in the case of the single crystals an 
insufficient amount of heat was transferred before reach- 
ing these peaks on the way up. The peaks were therefore 
not resolved in the up portion of the pulse. Lowering the 
maximum field of the pulse leads to a slight delay of the 
time when the peaks are reached on the way up, which 
improves the chance of resolving these peaks. A more 
drastic change occurs when a single crystal is crushed into 
powder. The resulting much larger surface/volume ratio 
improves the sample-to-bath heat flow substantially, and 
the peaks are also resolved in the up portion of the pulse. 

Another phenomenon seen in Fig. □ is that (for de- 
creasing H) a rapid rise of dM/dH occurs below about 
50 kOe. This rise is due to MST's from clusters in which 
the spins are coupled by DN exchange constants. The de- 
termination of these DN exchange constants is the main 
topic in the remainder of the paper. 
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FIG. 4: Magnetization traces for Zni_ I Mn I crystals with 
x=0.0056, 0.021 and 0.030 at T=0.65 K. These data were 
taken in dc magnetic fields. A (minor) correction for the 
lattice diamagnetism is included. 
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FIG. 5: Magnetization traces for x=0.029 measured at 
T=20 mK with Hlc (along the [1010] direction) and H || c. 
The magnetization was corrected for the lattice diamagnetism 
and normalized to its technical saturation value M 3 . 



DC MAGNETIZATION AT 0.65 K 



Figure 01 shows magnetization curves at 0.65 K, mea- 
sured in dc magnetic fields up to 170 kOe. As already 
noted, there was no hysteresis in any of the data taken 
in dc fields. The curves in Fig. 0] exhibit the expected 
behavior pS. Above 50 kOe the magnetization M shows an 
apparent saturation ( "technical saturation" ) . The appar- 
ent saturation value M s is lower than the true saturation 
value Mo. The latter is expected to be reached in fields 
substantially above 170 kOe. The relation between M 3 
and x was used in Sec. lIIII as one of the three methods of 
determining x. An expanded view of the upper portion 
of each of the curves in Fig. 0] does not show any MST 
between 50 and 170 kOe. (The magnetization change as- 
sociated with the small peak Q in Fig. [5] is estimated to 
be about 0.2%. This small change was not resolved in any 
of the dc data for x < 0.03). The absence of detectable 
MST's in the dc data between 50 and 170 kOe indicates 
that all MST's from DN pairs occur below 50 kOe. 



A feature of Fig. 0]which is most obvious for x=0.030 
is a magnetization "ramp" ending slightly above 40 kOe. 
Magnetization ramps are produced by the coalescence of 
broadened MST's^ The ramp ending just above 40 kOc 
is due to the coalescence of MST's from clusters involv- 
ing the largest DN exchange constant, defined earlier as 
J( 2 ). These MST's were not resolved at 0.65 K, but were 
resolved at 20 mK. 
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FIG. 6: Expanded views of (a) the low field portion, and (b) 
the high field portion of the magnetization traces in Fig. 



VI. DN EXCHANGE CONSTANTS AND 
SINGLE-ION ANISOTROPY FROM 0.02 K DATA 

A. Overall View of the dc Magnetization at 20 mK 

Figure [S] shows 20 mK data for :r=0.029, taken both 
with H || c and Hlc. These data are normalized to 
the technical saturation value, M s , and are corrected for 
the lattice diamagnetism. Expanded views of portions of 
these data are shown in Fig. |SJ The main features are: 

1. M rises quickly at low fields. This initial fast rise is 
typical. It is mainly due to the alignment of singles 
(Mn 2+ ions with no significant exchange coupling 
to other Mn 2+ ions)i£ 

2. The initial fast rise of M is followed by a ramp 
which ends near 5 or 6 kOe, depending on field 
direction. This ramp is shown more clearly in the 
expanded view of Fig.^a). 

3. A second ramp, smaller in height but spread over 
a larger field interval, follows the first ramp. The 
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second ramp ends near 40 kOe. Just above this 
field, M reaches technical saturation. 

4. The magnetization depends on field direction. This 
anisotropic behavior is more obvious at low H. 

Each of the two ramps is due to the coalescence of 
broadened MST'siS A well defined ramp usually corre- 
sponds to one series of MST's. The high-field end of 
such a ramp is near the last MST from this series, and it 
can be used to estimate the relevant exchange constant 
J. A more accurate value for J can be obtained if the 
MST's on the ramp are well resolved. 

The end of the second ramp, near 40 kOe, leads to 
the estimate jW/fegW- 0.5 K for the largest DN ex- 
change constant. The end of the first ramp, near 6 kOe, 
gives j' 3 '/fcB~— 0.08 K. To improve on these rough es- 
timates it is necessary to examine the MST's which give 
rise to each ramp, taking into account the relevant weak 
anisotropic interactions. 



relevant intra-pair exchange constant decreases. That is, 
the effects caused by the anisotropics are very small for 
the MST's originating from NN pairs of either class; small 
but easily detected for the MST's (ramp) from pairs in- 
volving j( 2 ); and very pronounced for the MST's (ramp) 
from pairs involving j' 3 '. 

The single-ion anisotropy can be described by the 
hamiltonian 

DSl + ^Si + S^ + St), 

where the z direction is along the c-axis. From 
Electronic Paramagnetic Resonance (EPR) at 77 K, 
D/kB=— 31 mK and a//cB=0.3 mKi^ Because the term 
involving a is relatively small, it was neglected in the 
present work. Values of the dd anisotropy constant for 
different classes of pairs are given in Table [I] They are 
based on a simple model in which the two spins in a pair 
are represented by two points separated by r n . 



B. Cluster Models and Anisotropies 

Cluster models play a key role in analysis of MST's in- 
volving DN exchange interactions £ Detailed information 
about cluster models and their statistical properties is 
given in Ref. U In the present work, only the largest two 
DN exchange constants, and J^ 3 \ were determined. 
The cluster models that were used in the data analysis 
depended on the field range of the data: 

1. For fields below 8 kOe, both J< 2 ) and are im- 
portant. These exchange constants correspond to 
two classes of DN's, but the identity of neither of 
these two classes was known at the beginning of 
the analysis. Therefore, the cluster models used 
for this field range included: a) the two classes of 
NN's, and b) any possible two classes of DN's se- 
lected from the six classes listed in Table [I] 

2. Because the magnetization above 10 kOe is hardly 
affected by J' 3 ', the cluster models for the field 
range from 10 to 50 kOe involved only the one 
class of a DN associated with j' 2 ', and the two NN 
classes associated with Ji and J{. All six possible 
choices for the class for the DN were tried. 

3. Analysis of the data between 50 kOe and 90 kOe 
was based primarily on a model which included 
only the two NN exchange constants. Simulations 
of the magnetization curves showed that in this 
field range the effects of the DN exchange constants 
J*- 2 -* and J*- 3 -* were very small. 

In addition to the exchange interactions, two 
anisotropies were also included in the analysis: single- ion 
anisotropy, and dipole-dipole (dd) anisotropy. The effects 
of these anisotropies on the MST's (or ramp) from pairs 
become more pronounced as the magnitude |J| of the 



C. Results at 20 mK and their Analysis 

1. Objectives 

The two objectives of the experiments at 20 mK were: 
1) to determine the values of J^ 2 \ and D, and 2) 

to identify the DN classes associated with j' 2 ) and j( 3 \ 
The analysis consisted of a number of steps taken in se- 
quence. In what follows, each step in this sequence is 
outlined, and the results are summarized. The main as- 
sumption in the analysis is that the Mn ions are randomly 
distributed over the cation sites. 



2. Procedures for Identifying DN classes 

The general procedure of associating DN constants of 
known magnitudes with different possible classes of DN's 
involves comparisons between experimental magnetiza- 
tion curves and computer simulations. 2,5 Separate simu- 
lations are carried out for all competing possibilities for 
the DN classes. The six DN classes listed in Table [I] lead 
to 6x5=30 possibilities for the two DN classes associated 
with j( 2 ) and j( 3 \ (Interchanging the order of the two 
classes of DN's leads to a new possibility.) 

In the present work these laborious simulations were 
postponed, because a preliminary identification of the 
DN classes was possible based on a simpler, and more 
physical, procedure. Two favorable circumstances per- 
mitted the simpler procedure: a large difference in the 
magnitudes of and j' 3 -*, and the availability of data 
for a sample with a very small x. 

For the lowest Mn concentration, :r=0.0056, the ramp 
which ends near 40 kOe is due primarily to j' 2 ' pairs. 
To a good approximation, the total magnetization rise 
AM' 2 ' associated with this ramp is therefore related to 
the fraction of Mn ions which are in such pairs. This 
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FIG. 7: Magnetization trace for £=0.0056, obtained at 
T=20 mK with H _L c. The ordinate m is the magnetiza- 
tion M normalized to the technical saturation value M a . This 
figure illustrates the procedure used to estimate the contribu- 
tion Am' 2 ' to m from all j' 2 ' clusters, and the contribution 
Am (3) to m from all J (3) clusters. 

fraction can be calculated for any possible choice of the 
DN class associated with j' 2 ), so that AM' 2 ' for any 
choice can be calculated and compared with experiment. 
Inclusion of J' 2 ) triplets, in addition to j' 2 ) pairs, in the 
calculation of AM' 2 ) makes the comparison even more 
reliable. 

Using the statistical tables in Ref. 0, the contribution 
of pairs and triplets to AM' 2 ) was calculated for each of 
the six possible choices of the DN associated with j' 2 ). 
Only one of these choices agreed with the experimen- 
tal value of AM' 2 ). Thus, a unique identification of the 
cluster class associated with j' 2 ) was achieved. A similar 
procedure was also used for a preliminary identification 
of the DN responsible for j' 3 ) , but in that case two pos- 
sible DN classes gave good agreement with experiment, 
so that a unique choice could not be made. The full nu- 
merical simulations, for all the samples, carried out at 
the end of the analysis, confirmed the preliminary iden- 
tifications of the DN classes on the basis of the data for 
x=0.0056. 



3. DN Class Responsible for J (2) 

The procedure of extracting the experimental ratio 
Am=AM /M s from the data, for the ramps associ- 
ated with j' 2 ) and j' 3 ), is illustrated in Fig. [7| For 
the higher-field ramp, associated with j' 2 ), it gives 
AM' 2 )/M S =1.3% for £=0.0056. The theoretical value 
of AM (2) /M S is obtained by multiplying (AM (2) /M ) 
by (Mo/M s ). The first ratio depends on the DN class n 
that corresponds to j' 2 ). For each possible n, this ra- 
tio was obtained from probability tables for pairs and 



triplets that involve only the exchange constant for class 
n Q • The probabilities are based on a cluster model that 
includes only the two classes of NN's and the DN class n. 
The second ratio, (Mo/M s ), was calculated from the so- 
called NN cluster model^ in which only the two classes 
of NN's are included. 

For small x the calculated ratio of AM' 2 ) /M s increases 
as the coordination number z n increases. This depen- 
dence on z n is the key for identifying the DN class that 
corresponds to j' 2 ). For £=0.0056 the calculated val- 
ues of AM' 2 )/M S are: 1.39% for the only DN class with 
z n =2; between 3.1% and 4.2% for the three possible 
classes with z n =6: and 6.1% for the two possible classes 
with z n =12. The experimental value 1.3% agrees only 
with z n =2, which corresponds to n=4 in Table [I] In 
the ideal wurtzite structure such a DN is reached from 
the "central" cation by moving a distance r4=c=^/8/3a 
along the c-axis (see Fig.^l. In Table[I]the exchange con- 
stant J(4) for rt=4 is designated as Jg. It is noteworthy 
that this largest DN exchange constant is not for the DN 
that is closest to the central cation. The closest DN is of 
class n=3, at a distance r 3 —^/2a. Thus, the magnitudes 
of the exchange constants do not decrease monotonically 
with distance. This non-monotonic dependence on dis- 
tance was predicted by some theories; 11 ! 12 and has been 
observed experimentally earlier 



4- Values of and D 

The values for Jg and D were obtained from analysis 
of the well resolved MST's on the ramp associated with 
j' 2 Wg" (see part (b) of Fig.©. These MST's stand out 
more clearly in the derivative dm/ dH of the normalized 
magnetization m = M/M s , shown in Fig. [SJ For three of 
the four samples in this figure the data are for both H _L c 
and H || c. The dependence of the fields at the MST's 
on field direction is caused by the anisotropy. Because 
the c-axis is an easy axis (D<0), the spins in the pairs 
are aligned more quickly for H || c than for H _L c. The 
faster alignment for H || c is more apparent in part (b) 
of Fig. El 

All curves in Fig. [8] show the last four MST's from 
Jg pairs. For £=0.0021 the first MST is also seen, but 
only as a "shoulder" on the fast drop of dm/dH at low 
fields. The fields at the last four MST's are all above 
15 kOe. Values of both Jg and D were obtained from an 
analysis of these fields, based on a pair hamiltonian which 
included the exchange interaction due to Jg, the uniaxial 
anisotropy governed by D, and the dd-interaction. The 
latter was calculated using the value in Table U for rt=4. 

In the first step of the analysis, an approximate value 
for Jg was obtained from the field at the last (fifth) 
MST, and the EPR value of D (Ref. GJ) was adopted 
as the initial value. Later, the predicted fields at the 
last four MST's, for both field directions, were calcu- 
lated for many sets of (Jg, D). The best match with the 
experimental values gave: Jg/fce=— 0.543±0.005 K, and 
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FIG. 8: Numerical derivatives dm/dH of m = M/M s , ob- 
tained from the experimental magnetization traces at 20 mK. 
The results are for two field directions: H || c and H || [1010]. 
The latter direction is designated as H _L c. The sample with 
x — 0.021 has been measured only with H || c. The traces 
have been displaced vertically relative to each other, but the 
gain is the same. 




D/kB=— 0.039±0.008 K. The presumably more accurate 
EPR value for D/k B is -0.031 K. The difference may 
be related to our use of the dd-interaction constant given 
in Table [I] As noted, this constant was obtained from a 
simple model in which the spins in a pair are represented 
by two points separated by r 4 . 

5. Possible DN Classes Associated with j' 3 ' 

Possible assignments of the DN class which corre- 
sponds to j' 3 ' were made on the basis of the magnitude 
Ato^ in Fig. [7] The experimental results give Am' 3 ' = 
AM' 3 '/AM S «7.3% for x=0.0056. This Am' 3 ' was at- 
tributed to the combined magnetization rise from pairs 
and triplets involving this exchange constant. The in- 
clusion of triplets in the theoretical calculation of Am' 3 ' 
was more important than in the calculation of Am' 2 ', be- 
cause the triplets/pairs population ratio was higher. The 
triplets/pairs ratio increases with z n . The lowest possi- 
ble coordination number, z n =2, is for the neighbor class 
n=4 associated with j' 2 ' and Am' 2 '. 

Theoretical values of Am' 3 ' were obtained from cluster 
models which included the three largest exchange con- 
stants (Ji, J{, J3) and any one of the DN classes with 
either z„=6 or z n =12. The calculated values are: ap- 
proximately 3.1% for all three DN classes with z n =6, and 
approximately 6.0% for the two DN classes with z n =l2. 
On this basis, z n is equal to 12, which leads to two pos- 
sible neighbor classes: n=6 or n= 7 (see Table |U). That 
is, j' 3 ' is either Jg ut or J' A . Another possibility (unlikely, 
but cannot be ruled out entirely) is that j' 3 ' is associ- 
ated with two different classes, both with z n =Q, that just 



FIG. 9: Low field portion of the dm/dH data with (a) H || c 
and (b) Hie. Different experimental traces are displaced 
vertically from each other, but the gain is the same. Also 
shown are simulations for the two types of pairs involving the 
two classes of neighbors with z n =12. The heavy dashed lines 
are for pairs with exchange constant J4 between the central 
cation and a neighbor of the class n=7. The dotted lines are 
for pairs with J3 Ut , involving a neighbor of the class n=6. 
The simulations are for the actual temperature T=20 mK, 
and they include the single ion anisotropy, involving D, and 
the dd-anisotropy. 

happen to have very nearly equal exchange constants. 

6. Value of J (3) 

Figure [5] displays dm/dH data for both H || c and 
H 1 c, in the field range relevant for the analysis of j' 3 '. 
Due to the smallness of j' 3 ' , the effects of the anisotropy 
are important, as can be judged from the strong depen- 
dence of the results on the direction of H. The dm/dH 
traces for H || c, in Fig. EJa), do not show the regular 
sequence of peaks observed in Fig. [S] for the MST's from 
J^-pairs. However, all traces in Fig. E3a) show a peak 
slightly above 2 kOe. The best resolution is for the sam- 
ple with a;=0.021. In this sample, two additional peaks, 
just below and just above 4 kOe are also resolved. 

The results in Fig. ED were compared with simulations 
based on a pair hamiltonian which included j' 3 ' and the 
two anisotropics. The simulations used D/k-Q=— 39 mK, 
as determined in Sec. I VI (Til and the dd-interaction con- 
stant taken from Table [I] Because the latter constant 
depends on the DN class, simulations were carried out 
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for the two possible neighbor classes with z n =12, i.e., 
n=6 with Jf ut , and n=7 with 

For H 1 c, pairs involving neighbors of either class 
have different orientations relative to H. The differ- 
ent orientations lead to different energies from the dd- 
anisotropy, and they give rise to MST's at slightly differ- 
ent fields. Therefore, in the simulations for the perpen- 
dicular field direction, H || [1010], pairs involving DN's of 
either of the two possible classes were divided into groups 
Pairs in different groups had different orientations rela- 
tive to H. There were three such groups for the DN class 
n=6, with J$ ut , and two groups for n=7, with J' A . The 
simulated magnetization curve was obtained by adding 
the results from all the groups. 

The simulations were carried out using the actual tem- 
perature, T=20 mK. Non-thermal broadening mecha- 
nism, such as local strains that give rise to a spread 
of the exchange and anisotropy interactions^ were ig- 
nored. In these simulations the values D/ks = — 39 mK, 
and the dd- interaction constant from Tabled were kept 
fixed. The value of was adjusted to obtain the best 
match with experiment for both field orientations. The 
final dm/dH simulations are shown in Figs. (a) and 
(b) as the dashed and dotted lines. Obviously, the re- 
sults of the simulations depend on field direction. For 
H || c, the overall structure of the experimental dm/dH 
traces is reproduced by the simulations. However, be- 
cause some sources of line broadening were neglected in 
the simulations, the detailed structure is better resolved 
in the simulations than in the experimental curves. For 
Hie the individual MST's are resolved in the simu- 
lations, even after the different orientations of the pairs 
relative to H are included. Experimentally, however, the 
individual MST's are not resolved for H 1 c. This differ- 
ence is attributed, again, to the neglect of some broad- 
ening mechanisms in the simulations. A crude way of ac- 
counting for the neglected broadening mechanisms is to 
replace the actual temperature T in the simulations by 
a higher effective temperature T c ff. The minimum T c g 
which leads to unresolved MST's for H _L c is 65 mK. 

Some features of the experimental data in Figs. Ha) 
and 9(b) are sensitive to the mag nitude of J (3 > . For H II c 
these features include the field at the most prominent 
peak, and the field at the peak associated with the last 
MST. These peaks stand out most clearly in the trace 
for x=0.021. For Hie the field at the rapid drop of 
dm/dH, which is at the end of the ramp associated with 
this series of MST's, is sensitive to the value of j( 3 \ The 
value of J^ 3 * 1 was determined from comparisons of these 
experimentally observed features with simulations that 
used different values of . Assuming that is Jf ut , 
the results gave j( 3 VfcB=-0.074±0.005 K. The alter- 
native, J (3) =4 gave -0.082±0.005 K. The first choice 
gives a slightly better agreement with the data, but in 
our view the evidence is insufficient for concluding that 
the DN class is definitely J£ ut . 



7. Simulations with both DN exchange constants 

Simulations of the magnetization curves, in fields up 
to 60 kOe, were carried out in order to confirm the pre- 
liminary identifications of the DN classes corresponding 
to j( 2 ' and j( 3 \ The simulations were for all samples, 
in contrast with the preliminary analysis that was car- 
ried out only for the sample with the lowest Mn concen- 
tration, £=0.0056. The simulations used cluster models 
which included the two NN exchange constants ( J\ and 
J[), J^ 2 ) and J' 3 ', all having the values quoted above. 
Each cluster model was based on a specific choice of the 
two DN classes associated with J^> and . As noted 
earlier, there are 30 possible such choices. The results of 
the simulations were not sensitive to any change of the 
DN class n provided that the coordination number z n did 
not change. 

Because all anisotropics were neglected in the simula- 
tions, the comparison was made with "isotropic" mag- 
netization curves obtained from the experimental data 
using the relation 

M isot = (My + 2M x )/3, (2) 

where My is for H || c, and Mj_ is for Hlc. The widths 
of the MST's exhibited by Mi sot are larger than the ther- 
mal width at 20 mK for several reasons: First, Mj sot is 
an average over different groups of pairs with MST's at 
slightly different fields. Second, DN exchange constants 
smaller than were neglected. Third, variations of 
the J's, and of the anisotropies, caused by local strains^ 
were not included. To match the widths of the MST's 
on the ramp which ends at 40 kOe, the simulations were 
carried out using an effective temperature T e g=100 mK 
instead of 20 mK. This change has no effect on the iden- 
tification of the DN's responsible for j' 2 ' and J^; the 
only effect is to smooth the curves. Simulations using the 
actual temperature, 20 mK, lead to the same identifica- 
tions of the DN classes. 

Figure EH shows simulations for x=0.029. The "ex- 
perimental" curve represents M; sot . Each simulation is 
labeled as n=(i,j, k, I). This designation means that the 
four exchange constants J\ , J[ , , that were used 
in the simulation, correspond, respectively, to the neigh- 
bor classes n=i,j, k, I in Table [I] All the simulations as- 
sumed that the first two neighbor classes (i, j) were those 
of the NN's, that is, = (1,2) or (2,1). (These two 

alternative choices are related to each other by an inter- 
change of the NN classes assigned to J\ and J[.) Both 
choices lead to nearly the same curves in this field range. 
The best agreement is with the simulation n=(l,2,4, 6) 
in which J^=J^ and J^=J^ ut . This corresponds to 
one of the two possibilities which were identified earlier. 
The other possibility n=(l, 2, 4, 7), not shown in Fig. 1101 
leads to a very similar curve, so that a definite unique 
choice of the DN class for j' 3 ' is not possible. 

The same conclusions concerning the neighbor classes 
associated with ,/( 2 ' and were reached from com- 
parisons of the data for the other samples (x — 0.0056, 
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FIG. 10: Comparison between the isotropic "experimental" 
magnetization Mi S ot=(Mij + 2Mx)/3, (see text) and numerical 
simulations based on the four exchange constants associated 
with the neighbor classes n=i,j, k, I. The simulations assume 
a random distribution of the Mn ions, and an effective tem- 
perature T e ff=100 mK. 



x = 0.021 and x — 0.030) with simulations. Thus, the 
earlier identifications of the DN classes are confirmed. 
The good agreement between the data and the simula- 
tions also lends support to a random Mn distribution in 
the studied samples, which is the main assumption in the 
simulations. 
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